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The  adsorption  of  a  fluid  of  hard  spheres  of  diameter  o  onto  a  planar 
surface  containing  a  triangular  lattice  of  adsorption  sites,  spacing  a,  was  studied 
using  a  model  which  is  equivalent  to  a  lattice  gas  with  n-body  interactions  that  are 
related  to  the  n-body  correlation  functions  of  the  fluid.  The  present  paper  extends 
our  previous  work,  which  included  only  pairwise  interactions,  to  also  include 
triplet  interactions.  We  discovered  a  simple  but  accurate  analytical 
approximation  to  the  triplet  correlation  function  for  three  spheres  in  mutual 
contact,  which,  when  combined  with  a  Muller-Hartmann  and  Zittarz 
approximation  of  the  critical  point  in  the  equivalent  lattice  gas,  yields  the  estimate 
pco^  =  0.8409  of  the  fluid  density  at  the  critical  point  of  the  condensation  transition 
which  occurs  at  the  fluid-crystal  interface.  This  estimate,  which  includes  the 


effects  of  both  pair  and  triplet  correlations  in  the  fluid,  is  significantly  higher  than 
the  value  pco^  ■  0.6678  obtained  if  only  the  effects  of  pair  correlations  are  included 
in  the  calculation. 
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I.  Introduction 

In  the  present  paper  we  calculate  the  effect  of  three-body  correlations  on  the 
critical  point  for  two-phase  coexistence  in  a  model  for  the  adsorption  of  a  dense 
fluid  onto  a  crystalline  surface.  The  model,1  described  in  Sec.  II,  is  equivalent  to  a 
two-dimensional  lattice  gas  with  n-body  interactions  which  are  related  to  the  n- 
body  correlations  in  the  fluid  at  the  surface. 

In  a  previous  paper2  we  used  the  model  to  calculate  the  two-phase 
coexistence  curve  for  a  fluid  of  hard  spheres  of  diameter  a  near  a  smooth,  hard 
wall  containing  a  triangular  lattice  of  sticky  sites  with  a  first-neighbor  spacing  a. 
We  used  the  Kirkwood  superposition  approximation3  to  estimate  the  n-body 
correlation  functions  of  the  fluid  at  the  wall  as  a  product  of  pair  correlation 
functions.  Using  the  Percus-Yevick  (PY)  virial  approximation  to  the  pair 
correlation  function,4*5  together  with  the  exact  two-phase  coexistence  curve  for  the 
equivalent  lattice  gas  with  pairwise  interactions,  we  calculated  that  the  critical 
point  of  the  two-phase  coexistence  surface  occurs  at  pco3  =  0.7092.  This 
represented  a  significant  improvement  over  the  value  PgO3  =  0.4611  obtained  using 
the  PY  pair  correlation  function  and  the  mean  field  approximation  to  the 
coexistence  surface.6 

In  Sec.  Ill  we  use  the  exact  coexistence  curve  for  the  lattice  gas  with 
pairwise  interactions7  and  the  accurate  Camahan-Starling  (CS)  approximation  to 
the  pair  correlation  function5*8  to  yield  PgO3  =  0.6678,  which  is  slightly  lower  than 
the  value  obtained  using  the  PY  pair  correlation  function. 

In  order  to  assess  the  accuracy  of  using  the  Kirkwood  superposition 
approximation  in  the  present  model,  in  Sec.  Ill  we  introduce  a  higher  order 
superposition  approximation  to  the  n-body  correlation  functions  which  includes 
not  only  pair  correlation  functions  but  also  triplet  correlation  functions  for  a 
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triangle  of  spheres  in  mutual  contact.  The  triplet  correlation  function  for  three 
spheres  in  rolling  contact  has  recently  been  calculated0  using  the  Percus-Yevick  3 
approximation  (PY3),  which  includes  triplet  correlations.10  We  discovered  that 
the  numerical  values  obtained  for  the  triplet  correlation  function  when  all  three 
spheres  are  in  mutual  contact  are  very  well  approximated  by  the  simple 
analytical  form 

&3  =  4  -  7th-  7t]2-  2t|3 

(gf)2  "  4(i  -  n)3 

Here  g°  and  g°  are  the  triplet  and  pair  correlation  functions  for  spheres  in 
contact,  and  q  =  (xc/6)p a3  is  the  packing  fraction. 

Within  the  above  approximation  for  the  n-body  correlation  functions,  the 
adsorption  model  is  equivalent  to  a  two-dimensional  lattice  gas  with  both  pairwise 
and  three-body  interactions  which  is  itself  equivalent  to  an  Ising  model  with  two- 
spin  and  three-spin  interactions.  The  interface  method  of  Miiller-Hartmann  and 
Zittarz,11  which  approximates  the  interface  free  energy  of  a  two-dimensional 
Ising  model  by  considering  only  a  restricted  class  of  possible  configurations,  has 
been  used  to  calculate  the  critical  point,  at  which  the  interface  free  energy 
vanishes,  for  an  Ising  model  on  the  triangular  lattice  with  both  two-spin  and 
three-spin  interactions.12  For  the  equivalent  lattice  gas,  this  estimate  of  the 
critical  point  together  with  Eq.  (1)  yielded  the  critical  density  pco3  =  0.8409,  which 
is  substantially  higher  than  the  value  obtained  using  the  Kirkwood  superposition 
approximation  and  the  lattice  gas  with  only  pairwise  interactions.  Triplet 
correlations  in  the  fluid  are  thus  important  and  can  significantly  affect  the  phase 
behavior  at  the  fluid-crystal  interface  in  the  present  model. 
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II.  The  Model 

We  consider  a  three-dimensional  model  for  adsorption1,2,6  which  consists  of 
a  fluid  of  N  hard  spheres  of  diameter  a  in  a  container  of  volume  V  with  a  hard, 
planar  wall  at  z  =  -  a/2.  A  triangular  lattice  of  sticky  adsorption  sites  with  first 
neighbor  separation  a  is  located  on  the  wall.  The  coordinate  system  is  chosen  so 
that  the  contact  plane  for  adsorption  of  a  sphere  on  the  wall  is  at  z  =  0,  and  an 
associated  triangular  lattice  A  is  located  at  the  contact  plane  such  that  if  the 
center  of  a  sphere  is  at  a  site  of  A,  then  the  sphere  is  in  contact  with  a  sticky  site  on 
the  wall. 

The  potential  Us(r^)  for  the  interaction  of  a  sphere  i  at  rj  with  the  sticky  sites 
is  given  as 

e-PUVi)  =  i  +  X  I  8(ri-Rs),  (2) 

Rs€A 

where  p  =  (kT)_1,  6  is  the  Dirac  delta  function,  and  X  is  a  stickiness  parameter. 
The  positions  of  sites  on  the  triangular  lattice  A  at  the  contact  plane  are 
represented  by  Rs. 

The  Hamiltonian  for  the  system  can  be  written  as 

H  =  H0  +  Hs ,  (3) 

where  H0  is  the  Hamiltonian  for  the  system  without  the  sticky  sites  (the  smooth 
wall  problem),  and 

N 

Hs  =  2  Ufrj) . 
a  i=l 


(4) 
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We  shall  approximate  U(Rlt...,Rn)  as 
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U(RX . Rn)=  X  W  +  XW3, 

nn  A 


(10) 


where  W  is  the  potential  of  mean  force  for  a  pair  of  spheres  in  contact  on  a  pair  of 
sites  on  A,  and  W3  +  3W  is  the  potential  of  mean  force  for  three  spheres  in  mutual 
contact  on  a  triangle  of  sites  on  A.  The  first  sum  in  Eq.  (10)  is  over  first  neighbor 
pairs  of  sites  in  {Rj,...,Rn),  and  the  second  sum  is  over  triangles  of  first  neighbor 
sites  in  (Rx,...,Rn). 

This  triplet  approximation  to  the  potential  of  mean  force  is  thus  equivalent 
to  the  superposition  approximation 


gR(Kt . r„)  =  n  g°  n 

nn  A 


(11) 


where 


g2°  =  e-3W, 

g3  /(gj)3  =  e“PW3  .  (12) 

Here,  g3  is  the  contact  pair  correlation  function  for  the  smooth  wall  problem  for  a 
pair  of  spheres  on  a  pair  of  first  neighbor  sites  on  A,  and  g3  is  the  contact  triplet 
correlation  function  for  three  spheres  on  an  equilateral  triangle  of  first  neighbor 
sites  on  A.  The  superposition  approximation  of  Eq.  (11)  differs  from  the  Kirkwood 
superposition  approximation3  used  in  earlier  papers1'2’6  which  did  not  include 
triplet  correlations  on  the  triangles  of  first  neighbor  sites. 

Substituting  Eq.  (9)  and  Eq.  (10)  into  Eq.  (8)  and  changing  from  a  sum  over 
the  positions  of  labelled  spheres  on  A  to  a  sum  over  lattice  sites  of  A  yields 
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Z/Z0  =  I  e  e  e  ,  (13) 

(ti> 

where  is  the  occupation  number  of  lattice  site  i,  and  where 

Xp}(0)  =  e^.  (14) 

Thus  Z/Z0  =  E,  the  grand  canonical  partition  function  for  a  lattice  gas  of  spheres 
on  A  which  have  a  chemical  potential  p,  a  nearest  neighbor  pairwise  interaction 
W,  and  a  three-body  interaction  W3  on  triangles  of  nearest  neighbor  sites. 

III.  Estimates  of  the  Minimum  Density  for  Two-Phase  Coexistence 

The  fraction  of  sites  of  A  which  are  occupied  by  spheres  is  given  as  1>2*6 

e-i*r&taa-  ll5) 

The  simplest  approximation  is  to  assume  the  lattice  gas  has  no  lateral 
interactions;  i.e.,  W  =  W3  =  0.  Eqs.  (13)  -  (15)  then  give  the  Langmuir  adsorption 
isotherm13 


e— (16) 

1+Xpj(0) 

Without  the  presence  of  lateral  interactions,  no  phase  transition  is  predicted  to 
occur  at  the  contact  plane. 

The  next  higher  approximation  is  to  assume  that  there  are  only  pairwise 
lateral  attractions  in  the  lattice  gas;  i.e.,  W  <  0  and  W3  =  0.  The  two-phase 
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coexistence  curve  and  the  adsorption  isotherms  can  be  calculated  numerically 
within  the  mean  field  approximation.  This  approximation  predicts  the  critical 
point  of  the  condensation  transition  to  be  given  as6 

-<P  W),  =  2/3 

„  (17) 

(gg)c  =  =  1.9477  . 


Using  the  PY  virial  approximation  to  the  pair  correlation  function4'6 


1  + q/2 
(1-T1)2  ’ 


(18) 


where  q  =  (ji/6)po3  is  the  packing  fraction,  yields  the  estimate  pco3  =  0.4611  of  the 
fluid  density  at  the  critical  point  of  the  condensation  transition  at  the  contact 
plane.6 

For  the  case  of  pair  interactions  in  two  dimensions,  the  two-phase 
coexistence  surface  is  known  exactly,2,7  end  highly  accurate  approximations  to 
the  adsorption  isotherms  can  be  constructed.13  The  critical  point  occurs  on  the 
exact  coexistence  curve  for  the  case  of  pairwise  interactions  on  the  triangular 
lattice  at2,14 


(e$>,  =  3.  (19) 

The  PY  pair  correlation  function  of  Eq.  (18)  then  yields  the  estimate  pco3  =  0.7092 
for  the  fluid  density  at  the  critical  point  of  the  phase  transition  at  the  contact 
plane.2  This  estimate  is  a  significant  improvement  over  the  value  of  the  critical 
density  obtained  using  the  mean  field  approximation  to  the  coexistence  surface. 

If  this  minimum  fluid  density  for  a  phase  transition  is  calculated  using  Eq. 
(19)  and  the  highly  accurate  CS  pair  correlation  function,5,8 


the  estimate  pco°  =  0.6678  is  obtained,  which  differs  only  slightly  from  the  value 
obtained  using  the  PY  pair  correlation  function. 

The  next  higher  approximation,  which  is  the  central  focus  of  the  present 
paper,  is  to  assume  both  W  and  W3  are  non-zero. 

For  an  Ising  model  on  the  triangular  lattice  with  the  Hamiltonian  (sj  =  ±1, 

J3>0) 

k  =  J  I  SjS:  +  J3  I  SjS:Sk  -  h  I  Sj ,  (21) 

<ij>  <»jk>  i 

the  critical  point  for  the  ferromagnetic  (J  <  0)  case  was  calculated12  using  the 
accurate  interface  method  of  Miiller-Hartmann  and  Zittarz11  to  be  given  as 
exp(-4^J)c  s  3,  so  long  as  the  three  spin  interaction  on  a  triangle  of  sites  was 
bounded  as  J3  <  3  I  Jl/2. 

This  Ising  model  is  equivalent  to  the  lattice  gas  of  Eq.  (13)  with  W  »  4J  -  8  J3 
and  W3  =  8  J3.  The  critical  point  for  condensation  in  this  lattice  gas  is  then  given 
as,  letting  g0  s  gg/(g§)2, 

(fio>c  s  3  (22) 

so  long  aa  W3  <  -3W/4,  which  is  equivalent  to  g0  >  (gf)^4- 

For  the  case  of  three  spheres  in  mutual  contact,  g0  has  been  calculated 
numerically  using  the  PY3  approximation10  by  Attard  and  Stell.9  The  values  they 
obtained  for  several  fluid  densities  are  given  in  Table  l.9  Using  the  CS  formula  for 
g£,  it  is  easy  to  see  that  g0  >  (g^) y4  is  satisfied  for  all  densities  studied.  Using  the 
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CS  formula  for  g3  and  the  PY3  values  for  g0,  it  is  also  clear  that  W3  >  0. 

For  the  case  of  three  spheres  in  mutual  contact,  Attard  has  proposed  the 
simple  approximation15 

g0  =  +  1)/2  (23) 

We  used  the  PY  and  CS  approximations  for  g9  together  with  Eq.  (23)  and 
discovered  the  approximate  relationship 


g0CS  =  (SoPY  +  goPY3)/2  • 


(24) 


Eq.  (24)  then  yields  the  simple  closed-form  approximation  g0PY3  =  go*,  where 


go  = 


4  -  7q  +  7r)2—  2ti3 
4(1  —  T])3 


(25) 


Values  of  g0*  calculated  using  Eq.  (25)  are  given  in  Table  1  and  compare  very  well 
with  the  PY3  values.9 

Using  Eq.  (22)  and  Eq.  (25),  we  obtain  the  estimate  pca3  =  0.8409  of  the 
minimum  fluid  density  necessary  for  a  phase  transition  to  occur  at  the  contact 
plane.  This  estimate  differs  substantially  from  that  given  by  the  lattice  gas  with 
only  pair  interactions  and  indicates  that  the  inclusion  of  triplet  correlations  for 
three  spheres  in  mutual  contact  is  important  for  the  study  of  adsorption 
phenomena  using  the  model  described  in  Sec.  11.  Since  the  inclusion  of  triplet 
correlations  in  the  model  resulted  in  a  repulsive  (W3  >  0)  three-body  interaction  in 
the  equivalent  lattice  gas,  it  is  understandable  that  the  resulting  estimate  of  the 
minimum  fluid  density  necessary  for  a  phase  transition  is  larger  than  that 
previously  predicted  using  only  the  effc  ts  of  pair  correlations. 
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The  effects  of  triplet  correlations  for  other  configurations  of  three  spheres 
on  the  triangular  lattice  can  be  neglected,  however,  for  the  calculations  of  Attard 
and  Stell9  of  the  triplet  correlation  function  for  the  general  case  of  three  spheres  in 
rolling  contact  can  be  used  to  show  that  the  three-body  interaction  energies  in  the 
equivalent  lattice  gas  which  result  from  these  other  configurations  are  negligibly 
small.  As  a  result,  the  value  pco3  =  0.8409,  which  includes  the  correlations  of 
three  spheres  in  mutual  contact,  should  be  a  good  estimate  of  the  location  of  the 
critical  point  of  the  phase  transition. 

Monte  Carlo  computer  simulations  have  been  performed  for  a  similar 
model,14  one  in  which  the  sticky  sites  are  replaced  by  small  sticky  regions  in  order 
to  facilitate  the  adsorption  of  a  significant  number  of  hard  spheres  in  a  reasonable 
number  of  Monte  Carlo  steps.  These  Monte  Carlo  results  indicated  that  three- 
body  correlations  may  be  quite  important  at  high  densities.  This  is  certainly 
consistent  with  the  results  of  the  present  theoretical  investigation. 
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Table  1.  Comparison  of  g0*  values  calculated  using  the 
closed  form  approximation  of  Eq.  (25)  to  the  accurate 
numerical  values  of  g0PY3  obtained  by  Attard  and  Stell.9 
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po3 

So* 

go™ 

error 

0.1 

1.0730 

1.0713 

0.2% 

0.2 

1.1641 

1.1610 

0.3% 

0.3 

1.2796 

1.2753 

0.3% 

0.4 

1.4282 

1.4225 

0.4% 

0.5 

1.6228 

1.6164 

0.4% 

0.6 

1.883 

1.878 

0.3% 

0.7 

2.239 

2.238 

0.04% 

0.8 

2.738 

2.751 

0.5% 

0.9 

3.46 

3.52 

1.7% 

